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Abstract 



We describe the null geometry of a multiple black hole event horizon 
in terms of a conformal rescaling of a flat space null hypersurface. For the 
prolate spheroidal case, we show that the method reproduces the pair-of-pants 
shaped horizon found in the numerical simulation of the head-on-collision of 
black holes. For the oblate case, it reproduces the initially toroidal event 
horizon found in the numerical simulation of collapse of a rotating cluster. 
The analytic nature of the approach makes further conclusions possible, such 
as a bearing on the hoop conjecture. From a time reversed point of view, the 
approach yields a description of the past event horizon of a fissioning white 
hole, which can be used as null data for the characteristic evolution of the 
exterior space-time. 
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I. INTRODUCTION 



Numerical simulations of axisymmetric space-times have enabled construction of the 
event horizon traced out by the evolution of dynamical black holes Jl|-§[ . Besides confirming 
behavior dictated by the general laws of black hole dynamics |7|||, the simulations supply 
further insight which was not anticipated from analytic theory. In the case of a head-on 
collision, they supply the details of how the holes form and merge || . In the case of rotat- 
ing collapse, they reveal how an initially toroidal structure is compatible with topological 
censorship ||10|| . In this paper, we present a 3-dimensional analytic description of the event 
horizon for a multiple black hole space-time which, in the axisymmetric case, reproduces the 
qualitative features of the above simulations. 

As a stand-alone item, a black hole horizon is a null hypersurface whose cross-sectional 
surface area monotonically increases and approaches a finite limit in the future. The number 
of black holes contained at a given time is not conventionally defined in terms of such a stand- 
alone picture but rather in terms of the number of disjoint sections given by the intersection 
of the horizon with a Cauchy hypersurface In the approach we present here, in the 

case of the head-on-collision the notion of "two holes" arises intrinsically from a preferred 
slicing of the horizon based upon an affine parameter along its generating null rays. 

For this purpose, we consider the geometry of a null hypersurface M whose surface area 
has a finite asymptotic limit. The intrinsic and extrinsic geometrical properties of a null 
hypersurface cannot be described in terms of the same "3+1" formalism used for a space- 
like hypersurface. In particular, the degenerate 3-metric of a null hypersurface does not 
define an intrinsic covariant derivative, in contrast to the case of a space-like hypersurface. 



Dautcourt [11] has presented an alternative formalism for inducing an intrinsic geometry 



on a null hypersurface from the embedding geometry. Geroch JT2[ used a similar approach 
to describe null infinity X as a 3-dimensional null hypersurface detached from the physical 
space-time M. In Sec. |J, we take an analogous approach to treat the null hypersurface 
traced out by a black hole as a stand-alone geometric object. 

Numerical investigations fl|-|6| of event horizons have used the Cauchy initial value prob- 
lem to evolve a black hole space-time throughout a domain of sufficient extent (both in space 
and time) so as to include an apparent horizon which has become approximately (but not 
exactly) stationary. This implies that the marginally trapped surfaces defining the apparent 
horizon have almost stopped growing in surface area. The event horizon is then located by 
examining the null geodesies which pass in the vicinity of the last quasi- stationary marginally 
trapped surface obtained in the evolution. The approach adopted in this paper is rather 
different. We have developed a method that constructs a null hypersurface whose surface 
area has finite asymptotic limit and whose intrinsic geometry satisfies all other requirements 
for a non singular event horizon. Thus, from all intrinsic criteria, this null hypersurface 
represents an event horizon 7i. The results of this paper depend only upon the intrinsic 
geometry of 7i and not upon the properties of the embedding space-time (which will be 
presented in subsequent papers). Even so, in order to show that the results of this paper are 
physically meaningful, it is necessary to discuss whether there exists a space-time in which 
H can be embedded. 

There is a formal construction of a vacuum space-time based upon the characteristic 
initial value problem posed on an intersecting pair of null hypersurfaces [p^ -|T5|. Here we 
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consider two null hypersurfaces intersecting in a topologically spherical surface Sq. Intrinsic 
geometrical data must be given on the two null hypersurfaces, as well as a quantity, called 



the twist P | 13 1 , which is analogous to an extrinsic curvature in the Cauchy problem. We 
can apply this formalism to the case where So is a cross-section of the event horizon at a late 
quasi- stationary time. Then, referring to Fig. |1], TC intersects the ingoing null hypersurface 
J + at Sq- Characteristic data given on the portion of 7i and J + to the past of So, leads 
to a formal solution of the vacuum equations in the domain of dependence D~ to the 
past of So- To the extent that So hes in the late time stationary region of the horizon, 
J + approximates future null infinity. So, for a horizon describing a binary black hole, 
the physically appropriate characteristic data on J + would describe the outgoing radiation 
emitted during the merger, although any data satisfying the constraints would formally lead 
to a consistent vacuum space-time. (As in the Cauchy problem, Einstein's equations imply 
constraints on the characteristic initial data but these reduce to propagation equations along 
the null geodesies so that it is simple to isolate the unconstrained free data.) 

There are theorems establishing the existence of solutions to this double null initial 



value problem |16|--|T8||. Although global issues remain unresolved, these theorems guarantee 
existence in some neighborhood of the initial hypersurfaces. The double null version of the 
characteristic initial value problem is equivalent to the world tube - null cone problem [T9 



in the case where the world tube is null. There exists a stable, accurate and efficient 
characteristic evolution code pO| , ^I| (the PITT code) which evolves this initial value problem 
and which could be applied to construct a numerical solution to the situation described in 
Fig. [I]. This provides the local existence of a space-time satisfying Einstein's equations in 
the finite difference approximation. Ideally, we would like to construct a numerical solution 
throughout the domain of dependence D~ but that is more problematical. The evolution 
algorithm requires the foliation of D~ by a one parameter family of null hypersurfaces J v . 
However, referring to Fig. [I], such a foliation becomes singular at Jm in the portion of the 
horizon corresponding to black hole merger. Thus one could only obtain the post-merger 
space-time by this procedure. This problem is due to a coordinate singularity, not a physical 
singularity, and from a mathematical point-of-view the space-time is extendable to earlier 
times; but just how much earlier cannot be answered by a characteristic evolution. Thus 
the analytic event horizon for a black hole collision, which we present in this paper, can be 
used as characteristic initial data to construct a vacuum space-time (either analytically or 
numerically) covering some domain preceding the merger. 

A major motivation for the present study stems from a potential indirect application 
of the results to a calculation of the gravitational waveform radiated by coalescing black 
holes using the PITT code. Waveforms from highly nonlinear, highly distorted single black 
holes have already been obtained with this code |20|j2T[1 . In these simulations ||20|| , the two 
null hypersurfaces were chosen to be (i) a portion of the past horizon (white hole) of a 
Schwarzschild space-time and (ii) an outgoing null hypersurface emanating from a slice of 
the past horizon to X + , which (in order to introduce distortion) contains incoming radiation. 
The code has equal capability of carrying out such simulations with the static Schwarzschild 
white hole replaced by a dynamic past horizon corresponding to a white hole which is 
initially stationary and later fissions into a pair of white holes. In a time-reversed scenario, 
the outgoing waveform from this white hole explosion corresponds to the incoming radiation 
incident from past null infinity on the merger of a pair of black holes. While this is not the 
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correct physical prescription of initial conditions for black hole coalescence, if the system 
were linear this incoming radiation could be used to draw inferences about the outgoing 
radiation. Such linearity is of course not expected. Yet any means of obtaining a handle 
on the merger waveform is of current importance. In addition, a solution of this problem 
would unambiguously yield the outgoing radiation from a white hole explosion, a system of 
at least academic interest. 

The time reversal of the intrinsic geometry of the event horizon of a black hole, results 
in the event horizon of a white hole, and vice-versa. The construction developed in this 
paper is naturally expressed in terms of a white hole that bifurcates in the future, and thus 
our method is presented for the case of a white hole. The reinterpretation of our results, in 
terms of a black hole horizon, is straightforward. Thus we pose our investigation in terms 
of a white hole horizon TC, which constitutes a portion of a null hypersurface Af, with the 
property that its surface area decreases into the future and has a a finite asymptotic limit in 
the past. The null rays of Af leave Ti at points where they meet other null rays of Af. Such 
endpoints can occur at a "crossover" point, where initially non-neighboring rays intersect, or 
at a caustic, where neighboring rays intersect. These properties follow from the fundamental 
theorems of black hole physics. ]7|,|| 

Note that a crossover point may also be a caustic, as in the case of a spherically symmetric 
light cone and also in the case of the prolate spheroidal light cone considered in Section 
In addition, a portion of Af can continue smoothly across a non-caustic crossover point, as 
in the example of the oblate spheroidal light cone considered in Section [IV|. However, 7i 
must end at such a crossover point. 

In this paper, we do not calculate the gravitational radiation emitted by the system but 
study only the internal dynamics of Ti. In Sec. we specialize to the axisymmetric case 
which facilitates an analytic treatment of the endpoints of 7i. We find that the bifurcation of 
Ti in a white hole explosion has the same "pair-of-pants" structure (in a time reversed sense) 
observed in the numerical simulation of a head-on collision of black holes. Further features 
emerge, such as the ultimate fate of an "eternal" pants leg, the details of toroidal black hole 



formation in the vacuum case and a bearing on a strict version of the hoop conjecture [22 



II. 3-DIMENSIONAL AFFINE NULL GEOMETRY 

It is useful to consider Af as one of the null hypersurfaces in the double null initial value 
problem for the vacuum Einstein equations. As first shown by Sachs O], the evolution of 
the double null problem requires as boundary data the intrinsic conformal geometry on A/", 
i.e. a metric 7^ expressed in terms of an affine parameter u, up to the conformal freedom 
7 a {, — > fi 2 7ab- In addition, the intrinsic conformal geometry must be specified on a null 
hypersurface meeting Af transversely at some cross-section So, as well as the intrinsic 2- 
geometry and certain extrinsic curvature quantities of So- Here we restrict our attention to 
the intrinsic properties of Af and So- As emphasized by Hay ward |14]JT5[| , it is important to 



consider specification of an affine parameter u on Af as part of the data. 

This intrinsic data obey the Sachs equations [^J governing the expansion or contraction 
of Af. As a consequence, the data determine a unique metric 7^ from the conformal equiva- 
lence class on Af. Here 7 a b satisfies the degeneracy condition 7 a bn fe = 0, where n b is tangent 
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to the generators of Af. We choose n a to have the affine normalization n a d a = d u . 

Our aim is to use a special choice of Af as a stand-alone model of a white hole horizon 7i. 
We require that Ti be complete in the past and that its surface area have a finite asymptotic 
limit as u — > — oo. H ends in the future at points where its generators intersect, either at a 
caustic or crossover point. 

The Sachs equations can be derived by projecting the relevant components of the 4- 
dimensional Einstein equation into 3-dimensional form. In doing so, we must deal with the 
degeneracy of j ab in representing the counterpart of 4-dimensional covariant derivatives as 
operators on Af. The formalism adopted here makes explicit use of the affine structure of 
Af. 

We begin with the four dimensional description of Af as a null hypersurface embedded 
in a vacuum space-time with metric g ab and covariant derivative V a . On Af, the affine 
tangent to the generators satisfies the geodesic equation n b V b n a = and the hypersurface 
orthogonality condition rf a \J h rf\ = We make no assumptions about the behavior of n a off 
Af. We project 4-dimensional tensor fields into Af using the operator 

Pi = % + nf (2-1) 

where l a is the unique outgoing null vector field on Af which is orthogonal to the affine cross- 
sections and satisfies l a n a = —1. We can set l a = — V a w, where u is any smooth extension 
of the affine parameter to a field in the neighborhood of Af. 

The projected metric 7^ = P^P b d g c d is the pullback of the 4-metric g ab to Af. When 
restricted to the 2-surfaces determined by the affine foliation, the projected contravariant 
metric 7°* = P^P\g cd is the (unique) inverse of the pullback of g ab . 

We introduce the shorthand notation _L T\ for the projection (to the tangent space of 
Af) of a tensor T%. Thus, _L n a = n a and _L l a — l a . In addition, the following useful formula 
hold on Af: 







JL C 2 n g ab 

The last two of these equations can be verified using the commutation relation 

0=±[C n ,±]. (2.9) 

Our purpose is to rewrite the projected curvature components $ a fe =_L n c n d R ca bd =-L 
n c (V c V a — V a V c )nt, in a form intrinsic to Af. By applying the above formulae, we obtain 

$ab = \cllab ~ \l Cd {C nlac )C nlbd . (2.10) 



-L n a 


(2.2) 


jl r 


(2.3) 


_L C 1 

- 1 - >~"n l a 


(2.4) 


_L C n n a 


(2.5) 


-L V [a n fe] 


(2.6) 


£"n~1ab 


(2.7) 


£>nlab- 


(2.8) 
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This further simplifies by setting 7^ = R 2 h a b and ^ ab = R~ 2 h ab where h ab 'C n h ab = 0. (This 
can be achieved by choosing R 2 as the determinant of the restriction of j ab to the surfaces 
of the affine foliation.) Then, in terms of the shear tensor E a & = £ n h ab , 

I 1 
&ab = -C n (R h ab ) + h ab RC n R — -R h c E ac E fed . (2-11) 
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The Sachs equations follow immediately from Eq. fl2.11| ). Taking the trace of <& ab results 

n c n d R cd = - is 2 , (2.12) 

where R c d is the Ricci tensor and where S 2 = (l/2)h ab h cd 'E ac 'E b( i. Then in the vacuum case 
it follows that 

C 2 n R=-- A RT 2 . (2.13) 
The trace free part of Eq. ( p.llj ) yields 

V ab = ^£ n (R 2 Z ab ) - l -R 2 T?h ab . (2.14) 

where ty ab are projected components of Weyl curvature. 

We proceed further to decompose the shear in terms of its normalized eigenvectors p a and 
q a satisfying Y, ab (p h + iq h ) = T,h ab (p h -iq b ). Then h ab = PaPb + QaQb and E ab = T,(p a p b - q a q b ) , 
where p a + iq a = h ab (p b + iq b ) and satisfies p a C n p a + q a C n q a = and p a C n q a + q a C n p a = 0. 

The goal is to solve Eq. (|2.13 ) in a way consistent with a white hole horizon K on a 
portion of Af. Thus the Weyl curvature must be non-singular on TC. This requires that 
the contraction ty ab with a unit vector (normalized with respect to 7 a b) yield a non-singular 
scalar field on TC (including the endpoints of TC)] i.e. that 

q, = iTV + iq a )(p b + tq b )V ab = R- 2 C n {R 2 Y) + tZ(q a C nPa - p a C n q a ) (2.15) 

be non-singular. 



III. CONFORMALLY FLAT NULL GEOMETRY 

The Weyl smooth solutions of Eq. (|2.13| ) have a large degree of freedom corresponding 
to the outgoing radiation crossing M . In order to restrict this freedom we consider solutions 
whose null metric is conformal to that of a null hypersurface embedded in a flat Minkowski 
space-time. 

Consider then the flat space case, where we denote the corresponding fields on M as 
R, h ab , h a , it, etc. For convenience, we write F' = CnF for tensor fields F. Since ^> ab = 0, 
Eq. ( |2.14j ) implies (p a + iq a )(p a — iq a )' = and R 2 T, = a, where a 1 = 0. The conditions on 
the eigenvectors may be summarized by p' a = (H/2)p a and q' a = — (E/2)g . 

The focusing equation ( [2.13j) now integrates to give 
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R 2 = (Au + Bf - ia 2 , (3.1) 

where A' = B' = 0. We adjust the affine freedom in u so that R — > ft as u — > oo and so that 
the two caustics of R are placed symmetrically. Then Eq. (|3.1|) reduces to 

# = («+^)(fi-±*). (3.2) 

We choose conventions so that the a > 0, so that the caustic corresponding to the q principle 
direction is reached first, moving along a ray in the direction of increasing u. 
The dependence of the eigenvectors is determined to be 

where (P + iQ )' = 0. The resulting w-dependent family of 2-metrics, comprise the classic 
description of the 2-geometries generated by the parallel map |24| of a surface So embedded 



in Euclidean space. The parallel map consists of translations by the same distance Am along 
each normal to the surface (identical to Huyghen's construction for propagating a wavefront). 
Equivalently, by considering So to be embedded at time t = in a Minkowski space-time, 
the translation along each ingoing normal null direction through the time Am = t, generates 
a null-hyper surf ace foliated by constant time slices S t . 

In applying this construction, we choose S to be convex so that S t traces out the 
flat space wavefronts of a null hypersurface converging in the inward direction. From the 
point of view of the flat embedding, a is the distance between the two caustics generically 
encountered along each null ray. 

Given such a flat space null hypersurface, with the convexity property that its caustics are 
reached at finite u, we generate a curved space null cone with the same conformal structure, 
i-e. 7a6 = fi 2 7af>- We thus set R = QR and h a b = h a b- We do not require that the two affine 
structures agree and set n a = An a so that d u = Adu and E = AS. Our goal is to investigate 
the properties of the foliation St determined by translating So through the curved space 
affine time Am = t. 

The curved space focusing equation ( (2.13j ) now reduces to 

A'(Q'R + SIR!) + A(Sl"R + 2SY R!) = 0; (3.5) 
and the Weyl curvature, defined in Eq. ([2.14Q , reduces to 

q ab = °— (Sl 2 A)'(p a p b - q a q b ) (3.6) 
with the Weyl scalar, defined in Eq. (|2.15|) , given by 

. = (3.7) 

Sl 2 R 2 
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The goal is to solve these equations to construct a non-singular white hole horizon 7i. 
Then \1/ must be non-singular on 7i and the scalar fields Q and A must be smooth positive 
functions, except possibly at the endpoints of TC. In addition, the surface area of Tt must 
approach a finite limit as u — > — oo. 

We require that u and u approach — oo together at the same rate so that we may restrict 
the affine freedom in u by requiring that u' — > 1. The surface area function must have a 
finite limit R i?^ as u — > — oo, corresponding to an irreducible mass || = Roo/2. 
Then, from inspection of Eq. (|3.7|), Qu — > —Roo and A — > 1, as u — > — oo. We also assume 
that these conditions are uniformly satisfied along each null ray in terms of a 1/u expansion. 
This puts constraints on the fields A and Q which satisfy Eq. ( |3.5| ). In order to apply these 
conditions it is convenient to introduce the function F = AQ 2 . Then the smoothness of ^ 
requires that F' = at a caustic. Also, the asymptotic conditions on Q and A as u — > — oo 
require u 2 F — > R 2 ^ so that, u{logF)' — > —2. 

Equation (33) can now be rewritten in terms of F and Q as 

[/o 9 (nfl)]' 



We can generate solutions to Eq. ( p.8| ) by making an ansatz for f2 and then integrating 
to determine F. The above smoothness condition that F must satisfy at a caustic is then 
automatically satisfied if Q is smooth. The asymptotic conditions require that the ansatz 
satisfy OR -»• R^ + 0{l/u 2 ) and u 2 ({lR - 1) -> -(T 2 /24. 

In order for the resulting model to represent a non-singular white hole, additional con- 
ditions arise at a shear-free ray. Along such a ray, the focusing equation ( |2.13| ) implies that 
d^R = with solution R = C\ + C 2 u. Accordingly, R must be a constant along each ray 
since it approaches a finite limit R^ as u — ► — oo. Thus we must require that our ansatz 
reduce to f2 = R^/R along a shear-free ray. Such rays occur at any umbilical point of 
So where the two curvature eigenvalues are equal and a = 0. For surfaces of revolution 
the poles are always umbilical. Umbilics are a major factor in determining the qualitative 
behavior of the white hole model. Along a non-umbilical ray, the completeness of TC as a 
white hole model requires that the range of u extend to a crossover point or caustic, where 
it hits another ray and the hole terminates. However, along an umbilic, the white hole need 
not terminate and can extend to infinite u. As we shall illustrate, this is the mechanism 
which leads to multiple black holes in a w-foliation of 7i. 

The behavior at umbilics imposes further conditions on our ansatz. For example, suppose 
that the white hole terminates along a set of non-umbilical caustics u = —a/2 which has 
an endless umbilic ray on its boundary where a = 0. Then u must be finite along the 
non-umbilic set but approach oo as a — > 0. Also, vl = 1/A must have the same behavior 
since the caustic set, including its boundary, is reached at finite u. 

The simple ansatz 

2 

n = - Roo (u+ a y 1 (3.9) 

±2{p — u) 



satisfies all the above conditions if the parameter p is chosen so that p > a/yYS. Then 
Q > in the white hole region contained inside u < —a/2. Furthermore, integration of Eq. 
gives 
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(2M-5p-/i) 2 (WM+i) {6 - LU) 

and 

9 (2u-5p-p) 2 ^^ 

U (12u{u - p) - a 2 ) 2 (2u - 5p + /j) 2 Wf-i) ' ^ > 

where p = \/l3p 2 — a 2 is real and positive. The asymptotic expansion of the integral gives 

u = u -12 p\nu + C + 0(\), (3.12) 

u 

where C is the integration constant. 

On the caustic set u = —a/2, we have 

* = *°* ( ° + *P-l> \ **" (3 13) 

81 (<r + 2p) 6 (a + 3p) \a + 5p + pj ' V ' ; 

which is manifestly regular; and 

ul ^ + 2 P Yf* + S P + »Y\ (3 . 14) 
a \a + 5p-pj 

which displays the required singular behavior as a — > at an umbilical point and is otherwise 
regular. 

It is instructive to examine the behavior of the extrinsic curvature eigenvalues defined 
on H according to 

_ P a P b C nlab _ E d u R 

Kp ~ 2R 2 ~ 2 + IT {6Ab} 



and 

q a q b C nlah _£ aj? 
2i? 2 2 7? 



H <i *~>n lab 

S = 009. =-77 + — • ( 3 - 16 ) 



For our flat space model, k p = l/(u — a/2) and & 9 = l/(u + a/2). Then /C(p, g ) = A(re( Pi9 ) + 
£] -1 £}&fi). For the ansatz (|3.9|) , these reduce to 

_ a(12u{u - p) - a 2 ) (12(u - p) 2 + a(a - 4m + 2p)) (-2m + 5p - ^[OWd-i] 
Kp ~ 9(w-p)(2 M -a) (-2m + 5p + At )2[(2p/M)+i] ^ • ^ 

and 

-a(12u(u - p) - a 2 )(12(w - p) 2 + a(a + 4m - 2p)) (-2m + 5p - ^) 2 [(Wm)-1] 



9(m - p){2u + a) (-2m + 5p + u ) 2 [( 2 ^) +1 ] 



jr. (3.18) 



Both k p and K g approach as the white hole approaches "equilibrium" as u — > — oo. How- 
ever, for large negative it it is clear from examining the dominant terms in Eq's ( |3.17| ) and 
( P-18| ) that k p > and n q < 0. So, although the mean curvature of TC is negative in agree- 
ment with net focusing, Tt expands in the p principle direction. As illustrated in Sec. |VL the 
effect of this expansion as it becomes stronger near the g-caustic provides the characteristic 
shape of a bifurcating horizon. 
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IV. THE CONFORMALLY SPHEROIDAL CASE 



The preceding analysis describes the dependence of the null geometry of Ti on affine 
parameter along a ray. In order to obtain a model of an exploding white hole we examine 
the global dependence of the geometry on the angular coordinates parameterizing the rays 
of Tt. In the simplest case, S is a sphere. Then a = on all rays, h ab = q ab the unit 
sphere metric, R = —u and the ansatz Eq. (|3~9l) reduces to R = Roo so that the horizon is 
stationary. A conformally spherical white (black) hole is topless (bottomless) along all rays. 

Since all topologically spherical geometries are conformally related, the conformal null 
geometry of a stationary white (or black) hole is unique. What distinguishes a Kerr horizon 
from a Schwarzschild horizon is the initial geometry of Sq, which determines R^. In the 
Kerr case, R^ is the conformal factor relating the geometry of the Kerr horizon to the unit 
sphere. (The spin of the Kerr hole arises from an extrinsic geometric quantity which can 
also be freely specified on So). In the Schwarzschild case, R^ can chosen to be a constant 
(independent of angle). 

In the case in which Sq is a prolate spheroid (ellipsoid of revolution), we will show how 
the identical horizon structure arises, as found in the simulation of colliding black holes. 
Furthermore, in the oblate case, we will show how a temporarily toroidal horizon arises, as 
found in the simulation of rotating collapse. 

We start with the 2-dimensional surface S describing the spheroid 

^ + £=1, (4.D 
a z o z 

which can be alternatively described in coordinates y A = (9, </>), by the map 

x = a sin 9 cos (4-2) 
y = a sin 9 sin (4-3) 
z = bcos9, (4.4) 

with 9 G [0, 7r] and (f> G [0, 2tt). The intrinsic metric of Sq is 

g AB dy A dy B = (a 2 cos 2 9 + b 2 sin 2 9)d9 2 + a 2 sin 2 9d(f) 2 . (4.5) 

The determinant condition provides a way to define R 2 as 

R 2 = (4.6) 
det(q) v 1 

(with qAB the unit sphere metric in (9,(j)) coordinates). Therefore, R 2 = 
aVb 2 sin 2 9 + a 2 cos 2 9. A straightforward calculation provides the principal radii of cur- 
vature of iS>o, 

and 
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We consider Sq to be isometrically embedded at time t = in Minkowski space and 
identify the u foliation of 7i with the (ingoing) null hypersurface emanating from Sq. We 
invariantly identify the function a in our model of 7i as the difference between the principle 
curvatures of So: 

\b 2 - a 2 1 sin 2 9Va 2 cos 2 9 + b 2 sin 2 9 
a = \r e -r^\ = — . (4.9) 

We set u = uq on Sq. Then our convention that u — midway between the two caustics of 
TC allows us to invariantly identify 



(r e + r>) (2a 2 + (b 2 - a 2 ) sin 2 9)Va 2 cos 2 9 + b 2 sin 2 9 

This determines both the metric 7 a & and affine parameter u for the flat spheroidal null cone. 

The conformal factor for the metric ^ a b = ^ 2 7 a & of the curved space version is determined 
by the ansatz Eq.(|3^9l)_Here, in order to satisfy smoothness conditions, we require that the 
parameter p > a M J\/13, where a M is the maximum value of a attained on S . For a < b^/~3, 
the maximum occurs at the equator and % = \b 2 — a 2 \/a. For the highly oblate case with 
a > 6\/3, the maximum occurs between the equator and poles and aM = 2a 2 /(3V3b). 

To determine the curved space affine parameter, we set u = u Q on S Q . Along with the 
condition that du/du — > 1 as u — > — oo this fixes the remaining affine freedom in u, i.e. the 
ray dependent integration constant C in Eq. ( ^.12[ ). With this choice, u = u on Sq and 
du/du — >• 1 as u — >• — oo. 

Figure || illustrates the features of the flat space spheroidal null hypersurface Affiat, with 
horizontal lines corresponding to the foliation St given by u = Uq + 1. In the illustration, we 
suppress the rotational symmetry. As discussed below, when the spheroid is prolate (oblate), 
the crossover points X where Affiat pinches off is a space- like line (disc). The features of the 
conformally related curved space version M are quite similar when viewed with respect to 
the Minkowski foliation St- The chief difference is the effect of the conformal factor on the 
expansion and shape of the surfaces St, which produces a finite surface area as u — > — oo. 
The features of N with respect to the curved space affine foliation S t , given by u = uo + 1, 
are qualitatively similar to those for St at early times. The interesting black hole physics 
occurs near the crossover region of where the foliations St and St have topologically different 
properties. These are best illustrated by embedding techniques. 



A. Embedding 

Embedding diagrams constitute valuable tools for visualizing the intrinsic geometry of a 
curved 2-dimensional surface. By embedding the surface in a flat 3-dimensional Euclidean 
space, one obtains a surface with the same intrinsic geometry. The following technique was 
developed by Smarr, who applied it to the description of the Kerr black hole [p5fl . More 
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recently, it has been employed to analyze the event horizon of the head-on-collision of black 
holes [p6fl . Here we describe its application to our model. 

The first step is to introduce the angular coordinate rj = cos# which makes det{qAB) — 1- 
In (u, 7], 0) coordinates, the intrinsic metric of the horizon is 

lab dx a dx b = n 2 R 2 {f- l dr? + fdcf) 2 ) (4.11) 

where 

(i-M (412) 

This can then be transformed into (u,rj,(j)) coordinates by the substitution u = u(u,r],(f>), 
where u(u, rj, 0) is determined by integrating Eq. ( |3.11|) . The results of this paper are based 
upon carrying out the integral by means of a Taylor expansion in u about Uq up to 6th order. 
(No substantial change in our results were seen in going from 4th order to 6th). 

Now, one can isometrically embed this surface in a 3-dimensional Euclidean space with 
Cartesian coordinates x % by the map 

x 1 = F(r}) cos0, x 2 = F(t]) sin0, x 3 = G{rj), (4.13) 

where 



= y Q 2 R 2 / f - F 2 V . (4.14) 

The quantities F and G are used to display the surface in the familiar 3-dimensional flat 
space at a given instant of time determined by the u foliation. Moreover, one can monitor 
the embedding at different instants of time and produce an "embedding history" which 
shows the evolution of the surface's geometry. By suppressing the (f) direction one can stack 
<p = const cross-sections of the embedding of the St foliation in a three-dimensional fashion, 
with the vertical axis labeling t and the horizontal axes labeling F = \J (x 1 ) 2 + (x 2 ) 2 and 
G = x 3 . 



B. The pair-of-pants 

We first describe the prolate case b > a in which the crossover points X in the flat-space 
model are also a line of caustics with respect to the principle direction. (Thus the <fi 
direction corresponds to the q principle direction). In Fig. |2| the rotational symmetry has 
been factored out so a Minkowski time foliation S t of the underlying flat space spheroidal null 
hypersurface corresponds to horizontal lines. The effect of curvature focuses the conformally 
related null hypersurface Af and introduces an upward bulge in the S t foliation, which gets 
enhanced at later times to produce the slice S* at which the white hole bifurcates. The 
vertical time sequence in our figures corresponds to white holes but the figures can be 
turned upside-down to depict the corresponding scenario for black holes, in this case a black 
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hole merger. The points C represent the caustics at the poles which are reached at finite 
times in the flat model but correspond to infinite affine times in the curved model. 

Profiles of the embedding diagram of St at various stages are shown in Fig. [3[ Proceeding 
backward in time from the initial prolate spheroid Sq, the cross-sections form the sphere Si n f 
as t — > — oo. Proceeding forward in time, they form the surface S* (also indicated in Fig. 0) 
where the white hole is at the verge of fissioning. Note that the two white holes which are 
produced each have a sharp point at their inner pole. 

Fig. |] shows a time stacking of embedding diagrams of the St foliation to form an 
(inverted) pair-of-pants and Fig. ^ gives a cutaway view of the bifurcation. The main 
features of the pair-of pants agree with those found in numerical simulation of the head-on 
vacuum black hole collision ||, as described in Ref's |9|j26|]. However, the analytic nature of 
the present work allows us to draw the following further conclusions. First, for the vacuum 
case, the pair-of-pants is eternal along the two umbilical rays at the poles. However, the legs 
pinch off and shrink asymptotically since every other ray eventually reaches the crossover 
X at finite u. Also, referring to the discussion following Eq. ( |3.18|) , the principle curvature 
K q corresponding to the 9 direction is everywhere negative, except at the poles and in the 
limit u — > — oo where K q = 0. It is most negative along the equatorial rays, which gives rise 
to the bow-legged shape of the pair-of-pants. 



C. Toroidal horizons and the hoop conjecture 

In the oblate case a > b, the crossover points X in Fig. § (in which the orbits of the 
rotational symmetry have been factored out) correspond to the same crossover points as 
in the prolate case. The difference between the two cases is the orientation of the axis of 
rotation. Whereas X lies on the rotation axis in the prolate case (and thus determines a 
caustic line under revolution), in the oblate case X rotates to form a disc. Only the outer 
rim of the disc (generated by revolution of the equatorial points C in Fig. ||) consists of 
caustic points. 

Note that the induced metric jablx of the crossover disc is single- valued (except at 
the caustic rim C where it is singular); i.e, its value does not depend upon whether X is 
approached from the top or bottom. This is because (i) by construction of M as a null 
hypersurface embedded in Minkowski space, ^y a b\x is single- valued and (ii) the conformal 
factor Vt has reflection symmetry with respect to the equatorial plane. More generally, in 
the absence of symmetry, establishment of consistency conditions for a single-valued metric 
on the crossover surface would be more complicated. 

In the oblate case the 9 direction has the smallest radius of curvature (corresponding 
to the q principle direction). The umbilical rays at the poles cross before they caustic, so 
that the infinite umbilical stretch in the St foliation for the prolate case does not arise in 
the oblate case. Thus the event horizon is completed in finite affine time, in contrast to the 
prolate case. 

Another important difference also arises in this case. The umbilical stretch produces 
toroidal cross-sections of the horizon, rather than the two spherical sections arising in the 
prolate case. The details of the formation of the torus are best understood in terms of the 
affine displacement Au(9) between Sq and X, as a function of the ^-coordinate of the ray. 
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In terms of Minkowski time, the corresponding time displacement is [p~0| 



Then 



, , /sin 2 9 cos 2 9 , A 

^&{ff) = ^ S j — + — . (4.15) 



Au{9) = / u'du, (4.16) 

JUQ 



where u' is given by Eq. ( 3.1 1|) and Uq by Eq. ( 4.10 ). In the oblate case, Au has its minima 



at the equator but Au has minima at the poles and exhibits a monotonic growth towards the 
equator. As a result, the St foliation first touches X at the space-time point where the two 
polar rays intersect, creating a surface of revolution with the same double teardrop profile 
as S* in Fig. [J[ but now rotated about the vertical axis through the center. In successive 
cross-sections, St forms a torus (with sharp inner rim), which shrinks to a circle as the 
horizon terminates. The tidal deformation introduced in the St foliation of the curved space 
model is somewhat analogous to the pair-of-pants shape, except now the identifications of 
the suppressed rotational symmetry lead to a toroidal topology of S t for a period of time 
following the bifurcation. This regime corresponds to the scenario found in the numerical 
simulation of the collapse of a rotating cluster of particles |]|||l(|]. 

Euclidean embedding is not possible for this full sequence of toroidal white hole forma- 
tion. Similar results were previously noted by Smarr [^] in regard to the non-existence of a 
Euclidean embedding for high angular momentum Kerr black holes. The schematic profiles 
in Fig. | indicate the qualitative topological features of the evolution. 

The analytic nature of the present approach allows us to draw further conclusions. In 
the oblate case, the p principle curvature direction in which k p is positive corresponds to 
the direction, so that the equatorial circumference C is always larger than its asymptotic 
value 2ttR oq . This has important bearing on the hoop conjecture [22|| , which in its original 



formulation would require that C < 4nM. The exact nature of the mass M and of the 
inequality were purposely left vague in the statement of the conjecture for purposes of 
further mathematical refinement. If we identify M as the irreducible mass associated with 
the surface area AttR 2 ^ then M = and C > AttM at every finite u. The largest value 

of C occurs at the equatorial rim of the crossover disc where 



Cx = W ^ J farM] « 2.546[4vrM], (4.17) 
(V lo + 6) 

for the choice p = ou (which maximizes the result). Although our model should not be 
expected to provide the sharpest bound, this result suggests a significant restriction on any 
plausible version of the hoop conjecture. 



V. CONCLUSION 

We have shown that it is possible to treat multiple black or white holes via a stand- 
alone-model of the event horizon based upon constraint equations for the characteristic 
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initial value problem. In this paper, we have concentrated on the constraint governing 
the internal geometry of the horizon. Remarkably, this single equation produces such rich 
results. Even more interesting features should be expected for models conformal to flat 
space null hypersurfaces with more structure than the spheroidal case considered here. 

In subsequent work, we will extend the treatment to the constraint governing the extrinsic 
curvature. The boundary conditions provided by the solution of this constraint problem is 
the missing ingredient necessary to evolve the exterior space-time by means of an existing 
characteristic code. 
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FIGURES 




FIG. 1. A portion of the space-time prior to the coalescence of two black holes. The parameter 
v labels the advanced time on a family of incoming null hyper surf aces. D~ is the domain of 
dependence of characteristic data given on the event horizon H and on J + . 
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FIG. 2. Spheroidal (Mfi a t) and conformally spheroidal (N) null hypersurfaces: Factoring out 
the rotational symmetry allows the foliations to be depicted as lines from pole to pole. The 
Minkowski foliation indicated by Si is drawn horizontally. The curved affine foliation is indicated 
by S t . 
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FIG. 3. Embedding snapshots profiling fission of an initially spheroidal white hole So in terms 
of the coordinates F and G: The white holes are the surfaces formed by rotating the profiles about 
the horizontal G axis. 
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FIG. 4. The embedding history of the axisymmetric fission of a white hole into two white 
holes: The time slices St are horizontal and proceed upward as t increases.. The spatial origin 
of the embedding axes is offset from the center for clarity. The suppressed symmetry dimension 
corresponds to a rotation about the G axis. The history extends into both the future and past of 
the initial surface So, which lies approximately halfway up the picture. A time reversed view gives 
the pair-of-pants picture for the head on collision of two black holes. 
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FIG. 5. A downward view of the portion of the embedding history from Sq up to the verge of 
bifurcation where the light rays at opposite points on the equator are about to cross. 
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FIG. 6. Schematic profiles of an initially spheroidal white hole So going through a toroidal 
stage. The rotation axis is vertical. The heavy solid line depicts the shape just prior to formation 
of the torus. The inner rim of the torus (dashed line) has a non smooth edge. 
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